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The LABOCA survey of the Extended Chandra Deep Field 
South: Two modes of star formation in AGN hosts? 
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ABSTRACT 

We study the co-existence of star formation and AGN activity in Chandra 
X-ray selected AGN by analyzing stacked 870/im submm emission from a deep 
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and wide map of the Extended Chandra Deep Field South, obtained with the 
LABOCA instrument at the APEX telescope. The total X-ray sample of 895 
sources with median rcdshift z~l drawn from the combined (E)CDFS X-ray cat- 
alogs is detected at > 11(7 significance at a mean submm fiux of 0.49 ± 0.04mJy, 
corresponding to a typical star formation rate around 3OM0yr~^ for a T=35K, 
(3=1.5 greybody far-infrared spectral energy distribution. The good signal to 
noise ratio permits stacking analyses for major subgroups, splitting the sample 
by redshift, intrinsic luminosity, and AGN obscuration properties. We observe a 
trend of star formation rate increasing with redshift. An increase of star forma- 
tion rate with AGN luminosity is indicated at the highest L2-iokeV ^ lO^^erg s"-*^ 
luminosities only. Increasing trends with X-ray obscuration as expected in some 
AGN evolutionary scenarios are not observed for the bulk of the X-ray AGN 
sample but may be present for the highest intrinsic luminosity objects with 
L2-wkeV ^ lO^^ergs"^. This behaviour suggests a transition between two modes 
in the coexistence of AGN activity and star formation. For the bulk of the sam- 
ple, the X-ray luminosity and obscuration of the AGN are not intimately linked 
to the global star formation rate of their hosts. The hosts are likely massive and 
forming stars secularly, at rates similar to the pervasive star formation seen in 
massive galaxies without an AGN at similar redshifts. In these systems, star for- 
mation is not linked to a specific state of the AGN and the period of moderately 
luminous AGN activity may not highlight a major evolutionary transition of the 
galaxy. The change indicated towards more intense star formation, and a more 
pronounced increase in star formation rates between unobscured and obscured 
AGN reported in the literature at highest {L2-iokeV ^ lO^^'ergs"^) luminosities 
suggests that these luminous AGN follow an evolutionary path on which ob- 
scured AGN activity and intense star formation are linked, possibly via merging. 
Comparison to local hard X-ray selected AGN supports this interpretation. Star 
formation rates in the hosts of moderate luminosity AGN at z~l are an order 
of magnitude higher than at z~0, following the increase in the non-AGN mas- 
sive galaxy population. At high AGN luminosities, hosts on the evolutionary 
link/merger path emerge from this secular level of star formation. 

Subject headings: galaxies: active, galaxies: starburst, infrared: galaxies 
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Introduction 



Observing the co-evolution of active galactic nuclei (AGN) and their hosts is key to 
understanding the similar cosmic evolution of the space density of luminous AGN and of 
the star formation rate density. This co-evolution also has to lead to today's fossil relations 
between remnant supermassive black hole mass and the properties of the host spheroid. 
Directly measuring the star formation rate of a high redshift galaxy is however particularly 
difficult for AGN hosts, since the AGN proper will disturb and overwhelm the rest frame UV 
and optical spectral and photometric star formation tracers already at modest luminosities, 
unless the AGN is obscured. Furthermore, star formation in the host may be noticeably 
obscured in particular at high star formation rates approaching that of high-redshift ul- 
traluminous infrared and submillimeter galaxies. Infrared observations can thus play an 
important role for these studies. One approach that has been successfully applied to high 
redshift AGN is to use the mid-infrared polycyclic aromatic hydrocarbon ('PAH') emission 
features which can be separated fro m the strong AGN mid-IR continuu m emission by means 



of low resolution spectroscopy (e.g.. lHouck et al.ll2005l : iLutz et al.ll2005l ). However, even with 



the superb spectroscopic sensitivity of the Spitzer Space Telescope, this approach is limited 
for high redshifts to modest sample sizes. Alternatively, observations of the rest frame far- 
infrared/submm continuum have been used since the advent of the first sensitive submm 
photometers to measure star formation via the far-infrared/submm continuum from the as- 
sociated cool (T~35K) dust. This rests on the submm continuum being due to star formation 
in the host galaxy, with star formation dominating over the AGN heated dust emission at 
these wavelengths for all but the highest ratios of AGN luminosity to star formation in the 
host. Such a star formation dominance in the submm is possible because o f the steep de- 



cline of AGN dust emission towards far-infrared and submm wavelengths (e.g.. lPier fc Krolik 
I992I . and later torus models), while the SEDs of star forming galaxies have a pronounced 
far-infrared (FIR) peak. The assumption of star formation dominating the submm emis- 
sion is supported by several types of observations. For NGC 1068, the AGN in the local 
universe for which current spatial resolution is sufficent to separate the AGN from the host 
in the far- infrared /submi l limeter range, this assumption is directly supported by observa- 



tions (jPier fc Kroliklll993l : iPapadopoulos &: Seaquist 



1999l : lLe Floc'h et al.ll200l[ ). Note that 



NGC 1068 falls wel l into the range of AGN luminosities (intrinsic L2-iokeV ~ 10 erg s 



Colbert et al.l |2002| ) and far-infrared luminosities {Ljr ~ 10 L0) of the AGN found in 
deep X-ray surveys, such as the (E)CDFS AGN that are discussed below. Recent Spitzer 
spectroscopic studies using the polycyclic aromatic hydrocarbon (PAH) emission features 
in comparison to the far-infrared/submm emission lend further support to this assumption, 
indicating that, at Lagn/ Lfir ~ 10, the far-infrared emission of local Q SOs as well as of 



mm-bright high redshift QSOs is dominated by star formation in the host (jSchweitzer et al. 



-4- 



20061 : 
(e.g., 



Netzer et al. 


2007 


Evans et al. 


2001 



Lutz et al.ll2008l). plausibly fed by large reservoirs of molecular gas 



Solomon &: VandenBoutI l2005l ) . Here, we study a sample that is 
at or below this ratio of AGN and star forming luminosities (§3.2). Our sample does not 
reach the extreme values as some of the most luminous but mm-faint high-z QSOs, for which 
th e assignment of f ar-infrared emission to star formation is more uncertain (see discussion 



m 



Lutz et al.ll2008l ). Our sample is X-ray selected and thus not biased towards extremely 



obsc ured infrared object s as, e.g., some local ULIRGs or high redshift dust-obscured galaxies 



[e.g., iHouck et al.l 120051 ). for which the AGN contribution to the far- infrared emission may 
be significant. 

Measuring the submm emission and on its basis the star formation rate of AGN hosts 
thus plays an important role in constraining the evolution of AGNs. Such studies directly 
benefit from the improving depth and areal coverage of current submm surveys. Both, 
the study of local AGN and of the AGN/galaxy (co)evolution have lead to the suggestion 
of models in which intense star formation events and p owerful AGN activ i ty are physi- 



cally linked and sequentially occur in a s ingle object (e.g., ISanders et al.lll988l : lFabianlll999 



Granato et al.ll2004J : iHopkins et al.ll2006[ l. in a process closely linked to the hierarchical merg- 
ing of galaxies in the universe. In a nutshell, galaxy interaction followed by merging with 
associated gas inflow may trigger a powerful burst of star formation and subsequently feed 
the central black hole(s) of the merger, to produce a luminous AGN. The AGN may then 
quench the star formation event and shed the obscuring dust, and finally emerge as an opti- 
cally visible QSO. In this evolutionary picture, an obscured AGN would sample the earlier 
phases of the AGN activity and would be associated with more powerful star formation in 
the AGN host, and thus stronger submm emission. As already implied by the typically 
non- merger morphology of local moderate- luminosity AGN (Seyferts), such an evolutionary 
picture may not always apply and its applicability needs to be studied as a function of AGN 
luminosity and redshift. 

A different behaviour would be expected in the context of the successful unified picture 
of active galactic nuclei proper, in which the differences between obscured and unobscured 
types of AGN are the consequences of different viewing directions on intrinsically identically 
structured systems. In fact, for the AGN and its immediate dusty surroundings, for example 
in the form of a 'torus', emission would be expected to be equally strong or stronger in the 
face-on/unobscured directio n compared to the edge-on direction, eve n at long far-infrared 
or submm wavelengths (e.g.. |Pier fc Kroliklll992uNenkova et al.ll2008l ). However, because of 
the intrinsically weaker AGN/torus emission at far-infrared and submm wavelengths, such 
differences will be easily washed out by even a modest amount of star formation in the 
host galaxy. This would lead to the expectation of little dependency of submm emission to 
obscuration in the unified AGN picture. This 'unified' view is by no means contradictory to 
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the evolutionary picture. Considering both these perspectives only stresses that differences in 
submm emission found between different AGN types will to a large extent reflect evolutionary 
coupling of periods of AGN activity and host star formation, rather than AGN physics 
and orientation alone. If evolutionary signatures are found in the submm, the underlying 
mechanisms like merging have to be clarified from additional information like morphology 
or dynamics. Combining evolutionary and unified perspectives also emphasizes the need 
to separately test for possible evolutionary effects in populations of high redshift AGN of 
different luminosities, that may follow different paths of AGN and host evolution. 



surveys to the first generation of (sub)mm surveys 



20001 : iHornschemeier et al.ll2000l : iBarger et al.ll2001 



fe.g., 


Fabian e 


al. 


200d; 


Severenini et al. 


Almaini et al. 


2003; 


Waskett et al.ll2003) 



In general, no straightforward correspondence between typical sources frora these X-ray sur- 
veys and bright submm sources detected with SCUBA (jHoUand et al.lll999l ) was established, 
and the average flux from submm stacking experiments of X-ray AGN was found to be 
low e.g., SssoMm = 1-21 ± 0.27mJy for the Chandra Deep field North (CDFN) jBarger et al. 
200lh and S«.nj,^ = 0.48 ± 0.27mJy from the Canada-UK deep SCUBA survey (CUDSS) 



( IWaskett et al.l 120031 ) . with insufficient S/N of the stacks for an in-depth analysis as a 
function of AGN properties. Pointed submm foUowup of selected, e.g., hard/lurninous 



sources from the deep X -ray surveys resulted in a few detections (e.g., iMainieri et al.l 12005 
Rigopoulou et al.l l2009l ) but also some upper hmits. Conversely, SCUBA sources were 
often found to be asso ciated with X-ray faint AGN in the very deepest X-ray surveys 
( lAlexander et al.ll2005al Jbl). These faint AGN do not dominate the energetics of the SCUBA 
sources and were undetectable in earlier analyses that were not based on the ultradeep 
2Ms Chandra data. Their black hole masses appear rnodest compared to sirn ilarly massive 



galaxies and to more powerful AGN (IBorys et al.l 120051 : lalexander et al.l |2008|) . 



The study that has been perhaps most successful so far in establishing bright submm 
emission for an X-ray selected AGN population and finding trends with AGN properties 
used a different approach. Selecting luminous X-ray absorbed but optically bright AGN 
not from deep field Chandra or XMM X- ray survey data b ut from identification of X-ray 
brighter sources from the ROSAT survey, iPage et al.l (l200l[ ) were able to detect 4 of 8 X- 
ray abso rbed QSOs at Ss5n;/m, > 5mJy. Later comparisons with matched X-ray unabsorbed 
samples (jPage et al.ll2004j : IStevens et al.ll2005l ) provided evidence for a lower submm detec- 
tion rate of the X-ray unabsorbed objects. This difference in submm brightness between 
obscured and unobscured objects supports the evolutionary view , but quest i ons re main in 
particular when comparing the large submm fluxes for some of the lPage et al.l ( 1200 ll ) objects 
with the more modest success of submm follow up observations of luminous hard sources 
from deep fields. If submm emission follows the emission at other wavelengths, this might 



- 6 - 



simply reflect the brighter observed fluxes and larger AGN luminosities of the iPage et al. 
( I2OOII ) sources. However, it should be noted that these are selected as moderately X-ray 



obscured {Njj ~ 10 cm~ ) but optically unobscured (optical broad emission line = Type 1) 
AGN, unlike many of the heavily obscured AGN from deep X-ray fields, which show Type 
2 optical spectra l acking broad lines. Like the broad lines, their bright optical magnitudes 
( jPage et al.ll2001bl ) argue against a significant obscuration of the AGN in the rest frame op- 
tical/UV. Extending such studies to include more typical obscured AGN is clearly important 
and a main motivation of this paper. 

We here present a study of the submm properties of X-ray selected AGN in the (ex- 
tended) Chandra Deep Field South ( E)CDFS. Making use of a new submm map provided 
by LABOCA at the APEX facility (IWeiB et al.l I2OO9I ) as well as current X-ray data with 
substantially improved identification status and characterisation of the AGN, we can study 
the submm properties and hence host star formation rates as a function of AGN properties. 
Throughout the paper, we adopt an Qm = 0.3, = 0.7 and Hq = 70 kms^^Mpc"^ cos- 
mology. When distinguishing between moderate luminosity (Seyfert) and high luminosity 
(QSO) AGN we refer to an intrinsic luminosity limit of L2-iokeV ~ ^ "^^^ 
otherwise. 



10 ergs unless stated 



2. Results from submillimeter mapping of the Extended Chandra Deep Field 

South 



The excellent X-ray and multiwavelength coverage of the 0.1 square degree Chan- 
dra Deep Field South (CDFS) and the surrounding 0.3 square degree Extended Chandra 
Deep Field South (ECDFS ), in combination w ith the powerful mappi ng capabilities of th e 
LABOCA submm camera (jSiringo et al.ll2009[ ) at the APEX telescope flGiisten et al.ll2006l ). 
enables us to take a fresh look at the issues discussed in the introduction, making use of the 
improved observational resources. We use the LABOCA 870/im map obtaiii e d by the LESS 
(LABOCA ECDFS Submm Survey) consortium as described by lWeifi et al.l (|2009[ ). We use 
the beam-convolved v2.2 final map which includes a total of about 350 hours of observing 
time over an area of about 40' x 40' with an rms noise level of about 1.2mJybeam~^ in the 
in ner 30' x 30 ' . Thi s map and the catalog of 126 sources detected at > 3.7cr is presented 



m 



WeiB et al.l (120091 ). We also make use of the residual map obtained by subtracting these 



fWeifi et al. 


2009; 


Coppin et al. 


2005) 
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2.1. Samples of X-ray AGN 



To maximise the statistics and to allow us to draw meaningful conclusions for physically 
selected subgroups of X-ray selected AGN, we use X-ray based AGN samples from Chandra 
observations for both the ECDFS and the deeper but smaller CDFS, fo r a total of 8 95 X- 
ray sources. For the CDFS, we mainly use X-ray spectral propertie s of iTozzi et al.l (120061) 
which are related to and based on the original CDFS observations of Giacconi et al. 
Recently, the de eper 2 Msec Chandra data and catalog for the CD FS have become av ailable 
(ILuo et al.ll2008[ ) but physical modelling comparable to the level of iTozzi et al.l (120061 ) is not 
yet completed (Bauer et al., in preparation). We therefo re use the update d 2 Msec obser- 
vational data (o bserved flux e s, har dness ratios etc.), from lLuo et al.l (120081 ) where available 
but stick to the iTozzi et al.l (120061 ) results for physical properties derived from X-ray spec- 
tral fitting (intrinsic luminosities. X-ray o bscuring co l umn density Nh-, etc.). We have also 
added 94 CDFS sources that are new from lLuo et al.l (120081 ). Lacking X-ray spectral fitting, 
these were used only for the combined stack and when analysing properties as a function of 
redshift, no modelled properties are available for those. To reduce contamination by non- 
AGN sources we exclude here objects that are likely nearby normal/star forming galaxies. 
Discrimination between such objects and AGN is approximately possible by comparing their 
X-ray to their optical p roperties; specifically we have adopted a cutoff log^fxlfii) > — 1 



(e.g., iBauer et al.ll2004l ). for the ratio of the X-ray flux to the optical R-band flux. Here 
we use as X-ray flux the observed full band flux if available, othe rwise the larger of hard 
and soft band flux. For similar reasons, we exclude 14 sources from iTozzi et al.l (|2006[ ) with 



intrinsic rest frame Ivn.5- infcp.v < lO^^ergs" ^. Sources with such low X-ray luminosities will 
mostly not be AG N (e.g.. lBauer et al.ll2004j ). We have furthermore excl uded 20 sources f rom 
Tozzi et al.l ( 120061 ) that are not re-detected in the dee per 2Msec data of ILuo et al.l (|2008[ ). In 
particular, this cut includes the lGiacconi et al.l ( 12002| ) XID 618. While t his so urce, indicated 



in only one of two source extraction methods used by iGiacconi et al. 



with an interesting z=4.76 submillimeter galaxy (SMG) (ICoppin et al. 



an X-ray source is not confirmed by the deeper X-ray data of ILuo et al. 



in our analysis as a luminous obscured X-ray source as inferred by ITozzi et al. 



(2002), is coincident 



20091 ). its nature as 
. Inclu ding it 
J2006[ ) would 



(2008 



increase the differences between luminous unobscured and obscured AGN that we discuss 
in Section 2.2. The overall CDFS list has 396 X-ray sources out of which 302 have X-ray 
spectral fitting from lTozzi et al.l ( 120061 ). 



For ECDFS sources outside of the CDFS, we use the X-ray catalog of iLehmer et al. 



(I2OO5I ). supplemented with currently available identification and red shift informati o n (Si l- 
verman et al. 2009, in preparation), but not yet including results of iTreister et al.l ( l2009l ). 
Again, we have imposed a log{fx/ fB.) > — 1 cutoff to reduce contamination by non-AGN. 
We avoid double counting objects in the iLehmer et al.l ( l2005l ) ECDFS catalog that are also 
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CDFS so urces. Object s in th e iGiacconi et al.l ( 120021 ) CDFS main catalog are identified by 
a fla g in iLehnier et al.l (|2005[ ). For these we use the CDFS data and spectral fitting only. 
The ILuo et al. (|2008l) data detect more sources at the edge of the CDFS that are in the 
Lehmer et al. catalog but not fiagged as also detected in the CDFS. We avoid double 



counting those by eliminating them via a 5" radius match and using the ECDFS X-ray data. 

For the CDFS sources of ITozzI et aP feooef) we a dopt the positions provided in that 
paper. For the additional sources from lLuo et al. (2008) we use their optical positions where 
given and X-ray positions otherwise. In the ECDFS, we use optical positions from identi- 
fications by Mainieri et al. (in preparation) when available and X-ray positions elsewhere. 
Note that the 27" FWHM beam of the LABOCA map, obtained after convolution of the raw 
map with its own beamsize of 19.2" ( IWeifi et al.ll2009l ). is large compared to the Chandra 
X-ray positional uncertainities even for the least favourable case of modest S /N faint sources 
at large off-axis angles, reducing the need for optical identifications for the stacking process. 
Typically, Chan dra X-ray positions are accurate to well below an arcsec (e.g.. Fig. 5 of 



Luo et all 120081 ). 



We adopt as redshifts for all the CDFS sources of lTozzi et all J2006h_ the value s given 
in that paper. About half of these are spectroscopic redshifts from ISzokoly et al.l (|2004| ) 
and other references, while the rest are photometric redshifts which benefit from the excel- 
lent multiwavelength coverage of the CDFS. For the ECDFS area and the identifications 
of Mainieri et al., we have used in order of preference (1) secure (e.g., two or more emis- 
sion lines or clear spectral features, Ca H&K) spectroscopic redshifts from Silverman et al. 
(2009, in preparation), (2) other spectroscopic redshifts from the compilation of Rafferty 
et al. (2009, in preparation), (3) photometric redshifts derived by Rafferty et al. (2009, 
in preparation) from a comprehensive compilation of multi-band UV to Spitzer-IKAC pho- 
tometry using the ZEBRA photo-z code and finally in a remaining 28 cases (4) photometric 
redshifts from the COMBO- 17 survey which provides accurate photometric redshifts using 
photometry in 17 pass-bands from 350 to 930 nm. Here we have used the latest v ersion 
of the COMBO-17 CDFS catalogue, following a calibration update JWolf et al.ll2008h . We 
use this dataset with two limitations. We consider only COMBO-17 sources with R< 24 
(Vega): at these magnitudes the errors on the photometric redshift estimates are expected 
to be less than l^phot — ^spec|/(l + ^spec) ~ 0.06. The COMBO-17 data for galaxies fainter 
than R= 24 (Vega) are too shallow for accurate photo-z determination of AGN. Further, we 
limit the use of COMBO-17 photometric redshifts to z< 1.2 because at higher redshifts the 
COMBO-1 7 estimates beconi e increasingly inaccurate due to the lack of NIR coverage (see 
Sec. 4.6 of IWolf et al.l ( 12004| )). We have wai ved this las t cons traint for objects best fitted 
with a QSO templates ( MCdass = 'QSO' in IWolf et al.l tooi)) for which the photometric 



redshifts are accurate at least to z^ 4 (see Fig. 18 of IWolf et al.l (|2004( )). In total, we thus 
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have redshifts fo r about three quarters of the X-ray targets outside the CDFS. For the new 
Luo et al.l (|2008[ ) sources we use redshifts from Rafferty et al. (2009, in preparation) where 
available. 



The combined CDFS+ECDFS X-ray based catalog has 895 sources with a median X- 
ray flux of 3 x 10~^^erg cm~^ s~^ in the observed soft and 10~^^erg cm~^ s~^ in the hard 
band. 748 sources have spectroscopic or reliable photometric redshifts in the range up to 
z~5 (median z=1.17), the median observed X-ray luminosity of those sources is 10^'^ergs^^. 
In the following we ca ll this the 'combin ed' sample, while with 'CDFS' sample we designate 
the 302 sources from iTozzi et al.l (120061 ) with X-ray spectral fits, remaining after the cuts 
described above (median z=1.04). 



2.2. Stacking procedure 



We derive average submm fluxes for a given source population by extracting both the 
map flux (in Jybeam"^) and the map rms noise, using bilinear interpolation to the source 
position in the beam- convolved LABOCA flux map and noise map. We then use the inverse 
variance weighted average of the fluxes measured over the stack. 

Stacking procedures for data from large-beam deep submillimeter maps that are ap- 
proaching the confusion limit at these wavelengths have to consider two effects. First, blank 
background and thus the 'zero point' of a sky image is hard to identify. At the RMS noise o f 
~1.2mJy of our data, integral number counts approach lO^deg"^ (e.g.. ICoppin et al.ll2006l ). 
which is only a modest factor away from the LABOCA beam density (2 x lO'^de g "^ for an 
area of vr HWHM"^ of the convolved beam); see also the discussion in IWeiB et al.l ( l2009l ) on 
the contribution of confusion noise in the LABOCA map. This means there is effectively no 
clean background sky. Second, for the large submm beam an elevated signal at the position 
of a stacked source may originate from the wings of the beam of a nearby unrelated bright 
submm source. Chopped beam patterns would cause additional complications but do not 
apply to the scanned LABOCA data that were obtained in total power mode. 

We have addressed this situation with simple Monte Carlo simulations, randomly plac- 
ing ~ 10^ Gaussian beams over a 1000x1000 pixel blank image and assuming a beam width 
(in pixels) equal to the LABOCA beam width. Input fluxes were distributed according to a 
simplified integral number count distribution with a po wer law slope -2 (approximating mea- 
surements of, e.g.. lCoppin et al.ll2006l : IWeiB et al.ll2009[ ) and extending 3 orders of magnitude 
down from the brightest source in the map, deep into confusion. If the simulated image is 
offset to a mean value and stacking experiments are done using the positions of subsets 
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of the input list, mean fluxes from stacks agree with the mean of the input fluxes within 
an error estimated from the standard deviation of the image, divided by the square root 
of the stacked sample size. The noise level of current submm maps and the fact that they 
do not reach below the knee of the integral counts implies a dynamic range between noise 
and brightest source of typically about an order of magnitude only. This reduces the effects 
of individual bright outliers on such error estimates. With a proper image zero point, the 
effects of the difficulty to define the background and of confusion thus cancel, and stacked 
fluxes can be used dir ectly. In the LAB OCA data, instrument noise is still (just) dominant 
over confusion noise (IWeifi et al.ll2009l ). its presence is again consistent with this stacking 
procedure. Fig. [1] shows the pixel flux histograms for the part of the LABOCA flux map that 
is within a factor 1.5 from the minimu m RMS of 1.07ni Jybeam~^ and for the residual map 
after subtraction of detected sources (IWeiB et al.l l2009l ). In stacking experiments we sub- 
tracted the respective means (0.154 mJybeam"^ and 0.072 mJybeam"^ for this part of the 
flux and residual map, respectively). This also ensures that zero average flux is returned from 
stacking random positions in the map. In deriving the stacked fluxes, we also calculate an 
inverse variance weighted stacked submap which, for stacks with significant detections, can 
be used to verify that the stacked beam is centered and reproduces the original beam of the 
map. Spatial clustering of the stacked population can in principle significantly affect source 
counts and stacking results. The effect on the stacks can be very important for instruments 
at similar wavelengt hs with much larger beams than that of APEX, like Planck- HFI (e.g.. 



Negrello et al.l 120051 ). It is negligible for the LABOCA beam at this wavelength (IBavouzet 
20091 ). and hence not considered in our stacking. 



We base our discussion on stacking the flux map but also provide stacking results for 
the residual map for comparison. In a few cases with poor statistics we explicitly discuss 
the effects of significant detections for the stacked population. Given the noticeable overlap 
between the subm illimeter population and weak X-ray sources representing moderate lumi- 
nosity AGN (e.g.. [Alexander et al.ll2003l . l2005bl ). exclusion of individually detected submm 
sources would bias our results. Future detailed identification campaigns and high spatial res- 
olution submm foUowup of the LABOCA survey will allow supplementation of this statistical 
approach with one based on the confirmed nature of individual sources. 

We quote below the stacked fluxes for various samples of ECDFS AGN. In assign- 
i ng errors we con sider that, while the instrument noise is gaussian to good approximation 
(IWeiB et al.ll2009l ). the pixel histogram of the LABOCA map is somewhat non-gaussian due 
to the effect of individually detected, as well as fainter, sources (Fig. [T]). Positive detections 
in stacks are hence more likely than for purely gaussian noise. We provide in Tables [T] and 
[2] two error measures for the mean flux of each stack. For a sample with N sources, the er- 
ror (Jmap provides the standard deviation from comparing the mean fluxes of many (>1000) 
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samples of N sources each, drawn at random spatial positions in the map. The value of amap 
allows assessment of the significance of a stack's detection. The error (yguhsampie provides the 
standard deviation from comparing the mean fluxes of many subsamples of N sources each, 
drawn randomly from the fluxes measured at the positions of our combined sample of 895 
AGN. This error is larger than amap since it also includes the spread in the properties of the 
AGN population. For that reason, we use it when assessing the significance of differences 
between subsamples of our overall AGN sample. We have also compared (Jsuhsampie to error 
estimates from bootstrapping into each subsample proper and found the latter estimates 
broadly consistent, but with large fluctuations due to the sometimes small subsamples that 
are used below. 



2.3. Stacking results for different AGN samples 



The stacked flux for the 895 X-ray sources from the combined CDFS+ECDFS sample 
is detected clearly at S87o^m=0.49±0.044mJy (11. la). The stacked image is shown in Fig. [2J 
A 2-dimensional gaussian fit to the stacked beam is centered at 2.2" ± 1.7" from the expected 
position and has a gaussian fit FWHM of 27" x 34" (±3"). These parameters support the 
correctness of the map reduction and the stacking procedure, given the nominal convolved 
beam of 27". Stacking the residual map after removal of all > 3.7cr LABOCA point sources 
again provides a clear detection at ~70% of this flux (Fig. [21 Tabled]). Use of the residual 
map will lower the average flux of a stacked X-ray population by excluding members of 
that X-ray population that are individually detected SMGs. It can also lower the average 
flux if a subtracted point source is dominated by an unrelated object but includes a weaker 
submm flux that is originating from the blended X-ray source proper. Then, the flux at the 
position of the nearby X-ray source will be oversubtracted by removing a point source with 
the combined flux. This is related to one part of the 'boosting' effect on the fluxes of low 
S/N source detections fror n such a map. This is avoided for our residual map in a statistical 



sense because 'deboosted' (ICqppin et al. 



deriving the residue (see also IWeiB et al. 



200 



5I) fluxes for the detections have been used when 



20091) . 



For redshift z~l, close to the median redshift of (E)CDFS AGN, and adopting a T=35K, 
/9=1.5 greybody far-infrared continuum shape for the star formation powered part of the 
SED, the total infrared luminosity corresponding to our 0.49 mJy LABOCA detection is ~ 
2.6 X IO^^Lq. The inferred mean star formatio n rate is ~ 27 Mp^ y r~^, assuming star formation 



dominated subr nm emission, the conversion of iKennicuttI (119981 ) and then multiplying by 0.6 



to convert to a 



Chabrierl (120031 ) initial mass function. Note that for z> 0.5 and given the 



negative K correction for submm emission (e.g.. Fig. 4 of iBlain et al.ll2002l ). the considerable 
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uncertainty of this estimate is mostly in the adopted temperature of dust heated by star 
formation (Factor <2 in luminosity for a 5 K difference), rather than in the difference between 
individual source redshifts and the redshift z=l adopted for the conversion. The average deep 
X-ray field AGN is thus residing in a moderately actively star-forming object, its luminosity 
placing it in the category usually called luminous infrared galaxies (LIRGs). This result 
certainly averages over a range of far-infrared luminosities but the LABOCA detection also 
in the stacked residual map argues that it is not only due to a few luminous outliers. Better 
characterizing this spread will be a task for the Herschel Space Observatory, the ~10mJy far- 
infrared SED peak expected for the adopted greybody is well within its capabilities. Given 
the > 11a detection of the combined sample, stacks for subgroups can still be detectable 
at good significance. We use such substacks in the following to probe for trends with AGN 
properties. 

We start with a simple splitting of the combined sample by redshift (see Table [1] for 
results of this and subsequent splittings of the combined sample). More than 80% of the 
sample have redshifts and are split into about equal groups below and above redshift z=1.2. 
There appears to be a trend towards higher star formation rates at higher redshift. A 
difference in mean submm flux between z<1.2 and z>1.2 sources is found at the 3. Oct level. 
As for other comparisons of AGN subsamples, we have adopted here the (Jsubsampie errors 
from Tables [1] and |2l Comparing z<1.2 to z> 2 sources gives a similar difference but 
only at the 2. lex level, due to the smaller size of the latter group. The sources without 
redshift assignment on average have slightly higher submm flux than any of these groups, 
consistent with the notion that a signi ficant fraction of them are located at high re dshift 



and remain more difficult to identify (lAlexander et al.l 1200 ll : iMainieri et al.l l2005al ) . We 



can also scrutinize this trend with redshift via a Spearman rank correlation test of the 
746 individual submm fiuxes with redshift measurements, rather than looking at binned 
averages. The correlation coefficient is a modest 0.108, not surprising given the significant 
noise on each individual fiux measurement, but the probability of exceeding this coefficient 
in the null hypothesis of uncorrelated data is only 0.003. Interpreting trends of submm flux 
with redshift in terms of total infrared luminosity assumes their proportionality due to the 
negative K-correction. This assumption starts to fail at z<0.5. However, changes in the ratio 
of submm flux and inferred infrared luminosity are small at the median redshifts of the bins 
discussed here (Fig. [3]). Deviations at z< 0.5 act in the direction of strengthening the trend 
of luminosity with redshift compared to the trend for submm flux. We illustrate this by 
reporting for those samples in Table [1] with complete redshifts also stacked IR luminosities, 
obtained from the luminosities of flducial T=35K, /3=1.5 greybodies matched to redshift 
and submm flux of each source. These stacked luminosities (see also lower panel of Fig. [3]) 
thus capture the effect of the redshift distribution of the samples, and conflrm its impact 



- 13 - 



to be small. Results from analysing redshift subgroups for the CDFS sample only, with its 
more complete identification status (see Table |2] for this and other results for this sample), 
agree with results for the combined sample in showing higher flux at higher redshift but with 
< 2cr significance of the trends at the given sample sizes. Of course, the trends with redshift 
(see also Fig. [3]) do not immediately imply evolution, because of the luminosity vs. redshift 
selection effects of the underlying X-ray surveys which are effectively flux limited. 

To test this further, trends with X-ray luminosity need to be explored. Roughly half 
of those sources in the combined sample with redshifts are not used here, because of the 
lack of X-ray spectral modelling for many of the ECDFS-only sources. For them, only a 
simple observed X-ray luminosity can be computed from the distance and observed X-ray 
fluxes. No significant trends of average submm flux can be observed with this observed X- 
ray luminosity, but this could also be due to the effect of variations in the obscuring column 
density erasing trends with intrinsic AGN luminosity. This aspect is bet ter addressed using 



only the CDFS sample, for which the spectral fits of iTozzi et al.l (120061 ) assign an intrinsic 
rest frame hard X-ray (2-10 keV) luminosity for all sources including hard/obscured ones. 
We find no noticeable differences between groups or with respect to the total sample when 
dividing at L2-iokeV = lO^^ergs"^ into a more and a less luminous group with about half 
of the sample each (Table [2l Fig. H]). This changes at the most luminous end: Sources with 
L2-iokeV > lO^^ergs"^ have more than twice the average submm flux, and differ from the 
average of the CDFS sample at 2.2a and from the sources with L2-iokeV < lO'^^erg s~^ at 2. 60". 
We further discuss this stronger star formation around the most luminous X-ray AGN below. 
Again the behaviour implied by the analysis of the stacks can be tested via a rank correlation 
of individual submm fluxes with log(L2-iofcey)- Given the upturn in submm only at the high 
X-ray luminosity end (Fig. H]) it is not suprising that no significant correlation is seen over 
the full range (N=302, C.C. =0.059, Significance 0.30), but correlation between submm flux 
and luminosity is found above L2-iokeV > lO'^^'^erg s~^ (N=98, C.C. =0.26, Significance 0.01). 

One of the most interesting parameters that is potentially linked by evolution to the 
level of star formation is the X-ray obscuring column density. Simple tests that can be 
done on the combined sample show no trend: Binning directly by X-ray hardness ratio 
as well as by crude estimates of the obscuring co lumn density from the location in the 



redshift vs. hardness ratio diagram (e.g.. Fig. 8 of ISzokoly et al.l 12004 ) does not show any 



noteworthy changes. We hence focu s again on the CDF S sample using the modelled X-ray 



obscuring column densities Nh from iTozzi et al.l (120061 ). As Table [2] and Fig. [5] show, there 



is no significant variation with obscuring column density even for this well-characterized 
sample and comparing subsamples that are individually detected at 3-5a. We here split the 
sample at two different column densities. First, Nh = lO^^cm"^ is often used as the limit 



distinguishing unobscured AGN from the larger number of obscured AGN (e.g., iTozzi et al. 
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20061 ). As a second test we broke the sample at Nh = lO^^cm ^ which would help identifying 



changes for the highest column density objects. These two cuts, as well as an intermediate 
one at = 3 x lO^^cm"^ that gives roughly equally populated bins above and below the 
threshold, do not reveal any significant trends with X-ray obscuring column density. For 
the Nh = 3 x lO^^cm"^ cut which has the most equally distributed statistics, the ratio of 
stacked submm fluxes for sources more/less obscured than the cu t is 1.1 ± 0.5, well below 
high ratios like the factor 4.4 found in the IStevens et al.l (120051 ) comparison of obscured 
and unobscured very luminous QSOs. Using submm emission as star formation indicator, 
there hence seems to be no clear trend of host star formation with nuclear obscuration for 
the typical Lo m^pv ~ lO'^^ergs^^ AGN. Recent attem pts using radio contiri uum emission 
(IRovilos et al.l 120071 ) and the [OII]A3727 emission line ( ISilverman et al.ll2009l ) to trace star 
formation agree in not finding such trends. 



2.4. Results for the most luminous CDFS X-ray AGN 

In the previous section, we have found for the full population of (E)CDFS X-ray AGN 
a trend of submm flux and star formation rate with redshift, a change with intrinsic X-ray 
luminosity only at the highest L2-iokeV ^ lO^^ergs"^ AGN luminosities, the border that 
is conventionally adopted between 'AGN' and 'QSO', and no significant change with X-ray 
obscuring column density. We will discuss the implications below but first specifically repeat 
the check for possible trends with obscuration at the highest AGN luminosities. This is 
motivated by the large range of luminosities covered, combined with the possibility that 
evolutionary paths may significantly differ between luminous QSOs and lower luminosity 
AGN. 



We first restrict the stacks which are formed via the iTozzi et al.l (120061 ) obscuring col- 
umn densities to AGN with high intrinsic luminosity L2-iokeV > 3 x lO^^ergs"^. A col- 
umn density Nh = lO^^cm"^ divides the sample roughly into half, but we also explore the 
traditional Nh = lO^^cm"^ separation. In both cases the more obscured sources appear 
brighter in the submm (Table [2]), but the difference stays well below 2cr significance in ei- 
ther case. As a next step toward the most luminous and reliable sources we further restrict 
to L2-^()kpv > KJ '^'^erg s~^ and sources with spectroscopic redshifts (redshift quality >1 in 



Tozzi et al.l (120061 )). Again we find the more obscured sources to be submm brighter (Tabled 
Fig. E]) with again a larger difference in submm flux if separating at Nh = lO^^cm"^, but low 
significance. Such a change in mean star forming rate with obscuring column density occur- 
ing at roughly Nh = lO^^cm"^ might indicate that the small variations with column density 
for the full CDFS sample may be dominated by the contribution of luminous AGN. Con- 
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sistent with this interpretation, a test sphtting the less luminous L2-iokeV < 3 x lO^^ergs"^ 
AGN at Nh = lO^^cm"^ further reduces the difference in submm flux between the two bins, 
compared to the full sample (Table 

An independent way of distinguishing unobscured from obscu red AGN is the op tical 
spectral classification. We use the highest luminosity AGN from iTozzi et al.l (120061 ). for 
which either -Lo.5-2A;eV or L2-iokeV is at least lO^^ergs"^. Among th ose, we restr i ct our selves 
to sources with spectroscopic redshift (redshift quality fl ag >1 in iTozzi et al.l (120061 )) and 
group them into optical Type 1 (BLAGN in the scheme of Szokolv et al. ( 20041), 10 sour cejll) 



or Type 2 (HEX, LEX or ABS classifications in the scheme of Szokolv et al.l (2004j ). 17 



sourcecl). Here we included XID 35 that was not c lassified bv ISzokolv et al.l (120041 ) but for 



which an optical spectrum from the VVDS survey (iLe Fevre et al.ll2004l ) is lacking obvious 
broad lines. Consistent with the analysis where obscuration was defined via the X-ray 
obscuring column density, the obscured optical Type 2 AGN are again submm brighter 
(Table [2]) than the unobscured Type 1 AGN, but the significance of the difference is low at 
just la. 

For both methods, the difference between submm fluxes of luminous unobscured and 
obscured AGN is stronger when stacking the flux map compared to stacking the residual 
map, suggesting a contribution of individually detected submm galaxies. Inspection shows 
that the 'hard' group deflned by Nh > lO^^cm"^ or by optical Type 2 in both cases includes 
two X-ray sources with S/N in the submm map above 3 at their position, because they are 
close to submm sources that are detected at > 3.7a. Both have extracted submm fluxes of 
>6mJy. The 'soft' group in contr a st ha s no such sources. These two SMGs are near the 
X-ray sources with iGiacconi et aP (l2002h XIDs 51 (z=1.097) and 112 ( z=2.940). VV ^ e have 
investigated these associations on the LABOCA map, the VLA data of iMiller et al.l (120081 ) 
and the FIDEL 24yum image (M. Dickir ison et al., in prepara tion). For XID 112, the position 



of the optical identiflcation for which ISzokoly et al. 



( I2OO4J ) determined the redshift agrees 



within 1" with a weak S^,4GH!■. =4:8nJy radio sou rce and a 24/xm source, and is only 2.4" from 



SMG LESS J33152.0-275329 (IWeifi et al.ll2009r ). This supports the association despite the 
fact this is a relatively complex region of the LABOCA map with two more submm sources 
within less than 1 arcmin. For XID 52, again the optical position agrees with a (93/iJy) radio 
source and a 24/im source, but the offset to LESS J33217. 6-275230 is 10.8", large (> 2a) 
for identifications of such a ~5(j submm detection and the 27" convolved LABOCA beam 



^ Giacconi et al.l ^20021) XIDs 11, 22, 24, 42, 60, 62, 67, 68, 91, 206 



IGiacconi et al.l (|2002|) XIDs 18 27, 31 , 35, 45, 51, 57, 76, 112 153 156, 202, 253 (changed redshift from 
Szokolv et al.l (|2004h in lXozzi et al.l (|2006l) . see also lRoche etaD ([2000)), 263, 268, 547, 601 
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(e.g., equation B 22 of llvison et al.l (120071 ) and the 6" positional accuracy estimated for the 
LABOCA map in lWeiB et aLl ( 120091 )). We maintain this source in our stacks given that there 
is no nearby alternative radio identification for the SMG, and the region is again complex 
with the next SMG detection almost blended, but this association is clearly uncertain. 

In our comparison of obscured and unobscured L2-iokeV > lO^^ergs"^ AGN we have 
found brighter submm fluxes by a factor ~3 comparing X-ray column densities above to below 
Nh = lO^^cm"^ and by a factor ~2 comparing optical Type 2 to Type 1. This is intriguing 
but the significance of the differences is too low in either case to claim a detection from our 
sample. The next steps in sample and field size and/or errors on the individual star forming 
rates will be needed for robustly confirming whether the stronger star formati on reported 



at hi gher AGN luminosity in X-ray obscured vs. unobscured broad-line AGN ( iPage et al. 



200ll . and subsequent work) also holds in the regime of the brightest X-ray sources found in 
deep surveys like the (E)CDFS. 



2.5. Results on individual literature X-ray sources 



CXOCDFS J033229.9-275106 (also called CDFS-202: iNorman et al.l l2002h at redshift 
of z= 3.70 has been consi dered a prototype of a luminous radio-quiet X-ray selected type 2 
QSO. ISturm et al.l (l2006l ) have obtained a deep Spitzer mid-infiaied spectrum of this source, 
detecting AGN continuum but not the 6.2/im rest wavelength PAH emission feature that 
could be detectable at this depth if the source were also hosting an extremely luminous 
SMG-like starburst. In accordance with this result, the LABOCA map does not show a 
detection at the position of CDF-S 202 (Sgyo = -1.01 ± l.lOmJy). 

From deep SCUBA fo l low u p of four heavily X-ray absorbed and X-ray luminous AGN in 
the CDFS. lMainieri et al.l ( 120051 ) report a Ssso = 4.8±l.lmJy detection for the z=3.66 object 
CDFS-263. The LAB OCA map gives a flux of Sgro = 1.88 ± 1.16mJy for this source. While 
we cannot support the iMainieri et al.l (l2005l ) result by an independent significant detection, 
the two measurements are still consistent within 2cr given the errors, and the source is 



l ikely among the submm brighter part of the X-ray AGN population. iRigopoulou et al. 



( I2OO9I ) publish results from an extension of this project to 8 sources. They do not detect 
at 850/im CXOCDFS J033229.9-275106, consistent with the LABOCA and Spitzer results 
reported for this source above. Stacking their 8 targets in the LABOCA map we find 
SsTo = 1.55 ± 0.41mJy, consistent with results for L2-iokeV 

> lO^^ergs"^ Type 2 AGN 

reported above but somewhat lower than the Ssr n = 4.0 ± O.SmJy obtained from scaling the 



variance-weighted mean of the IRigopoulou et al.l ( 120091 ) SCUBA results to 870yum 
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Koekemoer et al.l (|2004[ ) report the detection in the GOODS-S region of 7 'extreme 



X-ray/optical ratio sources' (EXOs) characterized by robust Chandra X-ray detections but 
optical nondetections to extremely low limits. All of these fall on the LABOCA map, none 
of them is individually detected and the stacked flux is 0.07±0.43mJy. While this nonde- 
tection at submm wavelengths excludes an explanation of EXOs by obscured AGN that are 
coexistent with extreme star formation, it is compatible with other possible spectral energy 
distributions of AGN at moderate to extremely high redshifts. 



3. Discussion 



Our data put strong limits on possible trends in submm brightness with obscuration for 
AGN with moderate luminosities, but are consistent with such a trend for luminous AGN. 
This result can be compared to previous studies. The first comparisons of deep X-ray and 
submm su rveys mentioned i n the introduction had too limited statistics to reliably address 
this issue. iPage et al.l ( 1200 ll ) used a different approach of selecting very luminous hard X-ray 
sources from the ROSAT survey and obtained significant submm detections in 4 out of 8 
z=l-2.8 obscured A GN, the weighted mean submra flux for all eight is 8350 = 4.4 ± O.SmJy. 
Subsequent papers (jPage et al.ll2004l : IStevens et al.ll2005l ) supported this result by matching 
to unobscured AGN of similar luminosity and redshift, and by extension of sample size, 
overall finding a ratio 4.4 in submm flux between X-ray unobscured and obscured AGN, and 
with the difference between submm detection rates of individual objects in the two groups 
significant at the 3-4(T level. 

Our LABOCA results grouping luminous (La-iofceV > lO'^^ergs"^) CDFS X-ray AGN 
by obscuration provide a ratio 2-3 in submm flux between obscured and unobscured objects, 
using separations by either Nh = lO^^cm"^ or by optical spectral type. However, due to the 
low significanc e we c an not rule out neither absence of a difference nor trends as reported 
by iPage et al.l (I2OOII ) and subsequent papers. One should note that, considering the ~8% 
correction from 870/im to 850/im flux densities for /3=1.5 optically thin dust at z=l-4, 
our averag e submm fl uxes of obscur e d sou rces are a factor ~2-3 lower than reported by 
Page et al.l (I2OOII ) and IStevens et al.l (120051 ) for their samples. Assuming that the increase 
of submm flux with X-ray luminosity at highest luminosities (Fig. H]) continues beyond the 
range up to L2-\okeV ~ 8 x lO'^'^ergs"^ covered by our sample, this difference r aay b e due to 
the typically higher -L2-ioxey ~ lO^^ergs"^ luminosities of the IStevens et al.l ( l2005l ) sample 
(see also Fig. |6]). Observations with better statistics at L2-iokeV ^ lO'^^ergs"^ will be needed 
to firmly establish whether in this regime not only the observed increase in submm flux 
occurs, but also a gap starts to open between submm properties of unobscured and obscured 
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AGN. 



An important point to note is that the IStevens et al.l ( 120051 ) obscured QSOs are X-ray 
absorbed (though with typical A^^^^ just above lO^^cm"^ not Compton thick) but optical 
broad line region (Type 1) objects. Stevens et al. are comparing optical Type 1 QSOs 
having low X-ray absorbing column density with Type 1 but higher X-ray absorbing column 
density. In contrast, we were first comparing groups defined by X-ray column density only, 
and in the second test comparing optical Type 1 to optical Type 2. Splitting our very small 
sample of 10 X-ray l uminous BLAGN further by X-ray column density to fully reproduce the 
Stevens et al.l ( 120051 ) approach did not show significant diffe rences of subgroups in this group. 
While our result is consistent with the IStevens et al.l ( 120051 ) finding of increased submm fiux 
for obscured luminous AGN, the detailed role of different definitions of AGN obscuration 
remains to be investigated with better sensitivity and statistics. 



Sturm et al.l ( l2006l ) used Spitzer mid-infrared spectroscopy to detect or put hmits on 



star formation in a sample of eight optical type 2 QSOs selected from deep X-ray sur- 
veys. These objects span a range of intrinsic luminosities i^o.5-iofcey = 10^^'^ to lO^^ergs^^ 
and X-ray obscuring column densities = 10^^'^ to > lO^^cm"^, with all but one above 
lO^^cm"^. We discussed the consistent LABOCA results for the z~3.7 source CDFS-202 
above. For the major par t of their sample, six sources at redshifts 0.205 to 1.38 (median 
~0.5), ISturm et al.l (|2006[ ) report one detection and five limits on mid-infrared PAH emis- 
sion with inferred star forming luminosities of the order 2 x 10^°Lq. This is lower than the 
~ 1.1 X IO^^Lq obtained from the mean LABOCA fiux for the z<1.2 redshift bin, averaging 
over all X-ray lumi nosities (Fig. [ 3 ]), con verting via our adoted T=35K, 13=1.5 SED and as- 
suming z=0.6. The ISturm et al.l ( 120061 ) sample is small and may have missed star forming 
objects, but another factor likely contributing to the difference is the extrapolation from 
the PAH measurem ents to total s tar fo rming luminosity which is a function of interstellar 
medium conditions. ISturm et al.l ( l2006l ) adopted a scaling factor based on the star forming 
galaxy M82. A scaling factor more similar to the one for high radiation fie ld intensity envi- 
ronme nts, as suggested for star formation in hosts of local type 1 QSOs by 



Schweitzer et al. 



( 120061 ). would increase the inferred star formation rates/limits of ISturm et al.l (|2006[ ) by a 
factor ~4 and bring them close to the typical z<l star formation rate inferred from the 
LABOCA submm fiuxes. 

Another modestly sized z^O.5 Q S02 sample studied spectroscopically with Spitzer was 
presented by IZakamska et al.l (120081 ). They observed 12 type 2 QSOs, 10 of which were 
selected from a large area SDSS sample primarily based on their optical Type 2 spectra, 
large [OIII] luminosities and bright mid-infrared continua. Two sources of their sample of 
12 entered the sample from different far-infrared selected programs. Compared to the Sturm 
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et al. QS02s they exhibit larger AGN mid-IR luminosities and, on the basis of X-ray to 
[OIII] ratios, possibly higher X-ray obscuring column densities. Zakamska et al. report PAH 
detections in 6 of 12 objects and infer a typical star formation luminosity of ~ 5 x lO^^LQfrom 
the median of detections and limits. Both of these Spitzer spectroscopic studies are broadly 
consistent with the LIRG-like star forming luminosities of luminous type 2 AGN found in 
our LABOCA study. The overlap in luminosity, redshift and other properties is not large 
enough for a direct comparison of methods. A combination into a single analysis is not 
straightforward, given the different uncertainties that are involved in extrapolation to far- 
infrared luminosity from either the mid-infrared or submm side. 



3.1. The co-evolution of X-ray survey AGN with their hosts: Two paths? 



The differences expected in the merger evolution scenario between star formation around 
obscured and unobscured A GN may be present only at the high luminosity end of our sample 
and have been reported by IPage et al.l ( 1200 ll ) at yet higher luminosities. For the bulk of 
lower luminosity L2-iokeV ~ lO^^ergs"^ sources we do not observe any significant trend with 
obscuration. This is likely indicating a different evolutionary path for these AGN and their 
hosts. 

The host properties of high-redshift Type 1 and luminous Type 2 AGN are difficult to 
constrain via common optical/near-infrared techniques, as the AGN often outshines the host. 
For the more accessible intermediate luminosity and Type 2 part of the population, several 
studies have concluded that z<l X-ray AGN are hosted by massive galaxies spanning the 
region from around the top of the 'b lue cloud' via the 'green valley' to the 'red sequence' in a 



color - absolute magnitude diagram (INandra et al.ll2007l : [Silverman et al.ll2008l : iTreister et al. 



20091 ). These are luminous host galaxies, with few of them fainter than rest fr ame absolute 
magn itude Mb = —20.5 and typically brighter than Mb = —21 ( see also iBarger et al. 



2008 



2003 ) . Photometric stel l ar mass analyses for z~l X-ray AGN hos ts (lAlonso-Herrero et al. 



Bundv et al. 


2008; 


Lehmer et al. 


2008: 



- stellar mass relation for the general z=0.7-l population (e.g.. ICooper et al.ll2008l ) suggests 
log(M* ]>10.5M(7^ massive host s. It is plausible from these results as well as from the local ev- 
idence (IKauffmann et al.ll2003l ) to assume that also the bulk of the z~l unobscured/luminous 
X-ray AGN population resides in massive host galaxies. 

Several studies have recently used mainly mid-infrared star formation indicators to con- 
clude that typical massive galaxies at redshifts 0.7-2 are almost constantly forming stars at 
considerable star formation rates, with the stellar mass normalized 'specific star formation 
rate' (SSFR) increasing with redshift. This is concluded from the presence of a fairly tight 
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Elbaz et al. 


2007; 


Daddi et al. 


2007) 



The tightness of this relation suggests a high duty 
cycle of star formation. For massive z~ 2 galaxies, a strong role of secular evolutionary pro- 
cesses compared to individual brief merger events is independently suggested by dynamical 



studies of rest frarne opt i cal/UV selected high redshift galaxies (e.g.. iForster Schreiber et al. 
20061 : ICenzel et allboosi : Ishapiro et allboosh . 



This developing picture of z~l-2 galaxy evolution is naturally complemented by our re- 
sults for the bulk of the X-ray sample: Strong trends of host star formation rate with AGN 
obscuration are lacking because such star formation rates are pervasive in galaxies of similar 
mass and redshift. The hosts are evolving secularly and star formation is not linked to a 
specific state of the AGN. The typical ~30Mq yr~^ estimated above for AGN hosts assuming 
star formation dominated submm emission compares well with typical st ar formation rates 
in z~l log(M*)~lO.5M0 galaxies jNoeske et all I2OO7I : IPaddi et al1l2007[ l. A further com- 
parative interpretation is currently not warranted given the uncertainty in the stellar masses 
of the AGN hosts and the fact that the comparison of star formation rates is subject to dif- 
ferent extrapolation effects. Star formation rates in these studies are based on extrapolation 
from the optical/mid- infrared while we extrapolate from the submm. This mismatch can be 
partly remedied by a comparison to submm fluxes (measuring the rest frame far-infrared and 
indirectly star formatio n rates) that were e stimated for optically selected galaxies from the 
LESS LABOCA survey. iGreve et al.l ( 120091 ) find for a K-selected Kvega < 20 sample stacked 
870/im fluxes of 0.17±0.01mJy (z<1.4) and 0.47±0.03mJy (z>1.4), similar but slightly be- 
low our trend for the AGN hosts, and again with a positive trend with redshift. The increase 
in AGN host star formation rate with redshift (Fig. [3]) is combined with a relation SFR to 
AGN luminosity that is flat over a wide range of moderate AGN luminosities (Fig. H]), i.e. 
SFR does not depend on AGN luminosity. This is consistent with the moderate luminosity 
AGN hosts indeed following the increase of SSFR with redshift of the general galaxy pop- 
ulation. Firmly establishing this behaviour will require studies with better SFR sensitivity 
for individual objects, removing the need to average over l arge stacks and better breaking 
redshift-luminosity correlations in the parent X-ray sample. [Silverman et al.l ( 120091 ) use the 
AGN-subtracted [OII]A3727 emission line to compare star formation rates in 0.48 < z < 1.02 
AGN and inactive galaxies of same stellar mass at same redshift. For moderate luminosity 
42 < /o5f(Lo.5--iofcey) < 43.7 they find indistinguishable SFR distributions in full agreement 
with the submm result. When including larger AGN luminosities a SFR excess is indicated 
in their data. If robust to the increasing technical difficulties of measuring mass and [Oil] 
SFR for luminous AGN, this could indicate the onset of star formation enhancement due to 
e.g. merging. 



Our interpretation that most of the moderate luminosity X-ray survey AGN are hosted 
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by massive secularly evol ving galaxies is consistent with morpho logical analyses. i/S'T studies 
of z~l X-ray AGN hosts (iGrogin et al.ll2005l : iPierce et al.ll2007l ) find the hosts to typically be 
bulge-dominated galaxies with only a modest fraction of hosts showing clear morphological 
evidence for recent major mergers, such as strong asymmetries. 

Analysis of deep X-ray surveys has clearly shown a difference in the redshift evolution of 
high and lower luminosity AGN, with the comoving density of lower luminosity AGN peaking 
at lower redshift than the z~2 'Quasar epoch'. Since cosmic halo merger evolution can be 
roughly matched to the evolution of quasars but not to the evolution of lower luminosity 
AGN, this has been interpreted in terms of a difference between merger driv en accretion at 



high AGN luminosities and more secular evolution at lower luminosities (e.g.. iHasingerl 12008 
Hopkins &: Hernquistll2009l . and references therein). 



All these lines of evidence place the bulk of the AGN population detected in deep X- 
ray surveys on a relatively gentle and secular evolutionary path. In these sources. X-ray 
obscuration may vary through orientation of the immediate AGN environment in a classical 
unified/torus picture, and possibly also with an additional contribution by obscuration on 
larger host scales, but obscuration is not intimately linked to the global star formation rate 
of the host. Only for the L2-iokeV ^ lO'^^ergs^^ QSO-like AGN, obscuration may coincide 
with high star formation, consistent with a classical merger evolutionary path. 



3.2. Star formation and AGN accretion rates 



Over cosmic time, the growth and merger rates o f black holes and their hosts have to 
estab lish the local black hole - bulge mass relation (e.g.. lMarconi fc Huntll2003l : iHaring fc Rix 
20041 ). In this context, it is interesting to compare for high redshift AGN populations the 
current accretion rate onto the black hole and the host growth via star formation with a 
relation that would establish on average the local ratio of black hole and bulge mass. Our 
results permit steps in this direction for X-ray selected AGN, within the constraints of results 
that are based on averaging over sizeable samples. 

In Fig. [6] we show the location of the CDFS stacks, grouped by intrinsic 2-lOkeV X-ray 
luminosity as in Fig. HJ now in a diagram comparing star forming and AGN luminosity. To 
obtain the AGN luminosities we have converted from median 2-lOkeV X-ray luminosity for 
each luminosity bin, assumed to be AGN domi nated, to monochrom atic 5100A luminosity 
using the luminosity-de pendent apx re l ation of ISteffen et al.l ( 120061 ) as cast into units suit- 
able for our purpose by lMaiolino et al.l ( 120071 ) . We have then converted to AGN bolometric 
luminosity adopting L5o/=7i^Li, (5 lOOA) (e.g.. iNetzer fc TrakhtenbrotI 12007 ). The star form- 



- 22 - 



ing luminosities are based on the stacked greybody luminosities for an adopted T=35K 
(3=1.5 SED at the indivi dual redshift of each source. For comparison, we add samples of low 



and high redshift QSOs (Netzer et al 



20051 : IPriddev et al.ll2003 



2007 



Omont et al 



Lutz et al.ll2008l : IPage et al.ll2004j : IStevens et al. 



20031 ). The local universe ratio of black hole to bulge 



mass ratio of 0.14% ( iHaring fc Rixll2004j ) can be converted into a formal 'steady growth lu- 
minosity ratio' for star formation and AGN LsF/LAGN=4:.7{0.1/ri) that is shown in Fig. [6] 
for a black hole accretion efficiency t]=0.1. 

The points shown for the CDFS X-ray AGN reflect the stack averages, i.e., the individual 
objects may scatter noticeably towards higher and lower star formation rates. Nevertheless, 
their location around the 'continuous growth' line is compatible with the picture of secular 
evolution outlined in the previous section. AGN accretion rates or star formation rates may 
here fluctuate to some extent with time, moving individual sources around this line in Fig.|6]in 
left /right and up/down direction, respectively, but the location of the population overall does 
not require to place them in any special evolutionary state that would be deviating from a 
long term growth. Only the highest luminosity L2-iokeV > lO^^'ergs"^ CDFS AGN approach 
the location of the local and high redshift QSO in the comparison samples. Compared to the 
continuous growth ratio, these most luminous X-ray AGN as well as the QSOs are growing 
their black holes at a much faster rate, indicating that the matching star formation is most 
likely spread over longer timescales. Again, we flnd that only the most luminous CDFS AGN 
match the possibly merger related evolutionary pattern of QSOs. 

The combination in Fig E] of the CDFS AGN w i th the more lu r ninou s high redshif t AGN 
obs erved in the (sublmm by IPriddev et al.l ( 120031 ) , lOmont et al.l (120031 ) , IPage et al.l (l2004j ) 
and lStevens et al.l (120051 ) conflrms the behaviour discussed in Sect. 2.1. Star formation in the 
hosts of modest luminosity high-z AGN seems to depend little on the exact AGN luminosity 

- reflecting the 'secular path'. These external samples, however, extend the increase in 
star formation above L2-iokeV ~ lO^'^erg s~^{Lagn ~ 3 x lO^'^erg s~^) that was already 
indicated in the CDFS data - reflecting the connection between star formation and AGN 
on the 'evolutionary connection/merger path'. The combined results for these samples and 
the CDFS AGN can be approximated by a simple relation Lsf = 10"^^-^^ + 10^^-'^ x L^^;^ 
(dashed line in Fig. EJ. Here, the smaller slope for the power law at high AGN luminosities 
compared to the slope 1 implied by the 'steady growth ratio' plausibly reflects that at the 
highest AGN luminosities the AGN growth (at the time of observation) is increasingly faster 
compared to the host growth. While this simple two component - constant plus power law 

- parametrization is a plausible reflection of the two growth modes, it should be considered 
illustrative and detailed parameters viewed with caution, given the sample selections and the 
fact that we here compare z~l AGN at moderate luminosities and z~2 AGN at the highest 
luminosities. 
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3.3. Comparison to the local AGN population 

The two evolutionary paths outlined above for high redshift AGN from deep X-ray fields 
imply a straightforward consistency check with the local AGN population. If the hosts of 
moderate luminosity AGN at z~l have star formation rates similar to nonactive massive 
galaxies at the same redshift, the star formation rates in the host of moderate luminosity 
AGN should follow the decrease of star formation on the general massive galaxy population 
towards redshift zero. 

To perform this check in a methodology as consistent with the LABOCA analysis 
of (E)CDFS X-ray AGN as possible, we have used a local unbiased 14-150keV extremely 
hard X-ray selected A GN sample from the 39 month Palermo Swift-BAT catalog (PSB, 



Cusumano et al.ll2009l see also Tueller et al. (2008) for an earlier BAT AGN catalog) in 
conjunction with the IRAS all-sky far- infrared survey. We selected (Shao et al., in prepara- 
tion) from the PSB survey sources classified as Seyferts, LINERs, quasars, and other AGN, 
explicitly omitting blazars. We excluded remaining objects with the possibility of a strong 
nonthermal contribution to the far-infrared on the basis of the NED SED, objects at galactic 
latitude \b\ < 15 and objects at redshift z>0.3 for which the IRAS 60/im band no longer 
probes the rest-frame far-infrared. For the remaining 293 AGN we used the IRAS Faint 
Source Catalog eO^m detections where available, otherwise we used Scaup! to obtaiu 60„m 
measurements for faint or individually nondetected objects. We calculated rest frame 2- 
lOkeV luminosities extrapolating from the BAT fluxes and the redshift, assuming an AGN 
photon index of 1.8, and infrared luminosities z/L^(60yum) in the observed frame. Because 
about 20% of the sample are individually undetected at 60yum and for consistency with the 
LABOCA stacking, we stacked the luminosities in seven 2-lOkeV luminosity bins spanning 
the PSB luminosity range with sufficient statistics in each bin. 

Star forming and AGN luminosities for these local AGN stacks are shown in Fig. [6] 
in direct comparison to the (E)CDFS results. At moderate AGN luminosities {Lagn < 
3 X lO^^erg s~^), the star forming luminosities of local hosts are about an order of magnitude 
lower than in z~l (E)CDFS AGN hosts. On the other hand, these local AGN appear to 
follow down to lower AGN luminosities the diagonal correlation of star forming and AGN 
luminosity that is reflecting the 'merger path'. T his is in full agree r nent w ith previous local 



work, for example the infrared-based results of iRowan-RobinsonI (119951 ) and iNetzer et al. 



(120071) coveri ng the regime of local QSOs, and modelling based on SDSS spectroscopy by 
Netzerl (120091 ). The latter traces a correlation down to even lower AGN luminosities < lO^^erg 



s ^ where star formation rates in our BAT stacks flatten out. We speculate this might relate 



'^http://scanpi.ipac.caltech.edu:9000/applications/Scanpi/index.htnil 
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to a more difficult disentanglement of weak AGN and moderately star forming galaxies 
galaxies in optical spectra than in very hard-X vs. far-infrared. 



MuUaney et al.l ( 120101 ) use Spitzer 70yum data to trace the evolution of the far- infrared 



to X-ray luminosity ratio of AGN to redshift z~2. They find no significant change of IR lu- 
minosity with redshift for luminous L2~iokeV = 10^'^~'^^erg s~^ X-ray AGN but an increase for 
more modest L2-iokeV = 10^^~^^erg s~^ AGN. While an observed wavelength of 70/im makes 
for a difficult AGN/host diagnostic at the high-z/high-1 end of that range, where AGN heated 
dust will enter strongly the observed band, these results agree with our submm-based finding 
that an increase of host SFR from local towards z~l occurs for modest luminosity AGN only. 
We note that the FIR luminosity increase by about an order of magnitude to z~l would 
be difficult to reconcil e with the alternative explanation of an AGN covering factor increase 
jMuUanev etldlboioh. since cove ring factors <0.1 would be required for local AGN. Another 
recent study (ITrichas et al.ll2009l ) addresses the starburst/AGN connection in luminous z~l 
SWIRE AGN. Being based on rare 70/im detections only it is not directly comparable to our 
work but consistent with the presence of elevated star formation in L2-iokeV ~ 10 
z~ AGN. 
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Our finding of elevated star formation in modest luminosity high z AGN is consistent 
with the change of typ ical star fornaatioii rates from z~l t o z=0 in the non- active mas- 
sive galaxy population f iNoeske et al]l2007l : Isibaz et al. l l2007l : lOaddi et all 120071 ). In Fig. il 
both local and z~l AGN are consistent with a 'secular' path with star forming luminos- 
ity independent of AGN luminosity, and an 'evolutionary connection' path with increasing 
star forming luminosity at high AGN luminosity. The evolutionary path, possibly linked to 
merging, appears to emerge locally at lower AGN luminosities from the locally lower level of 
general star formation. 

We have previously discussed the location of local and high redshift X-ray selected AGN 
and optical QSOs in the Lsf vs. Lagn diagram of Fig. |6l In the merger evolutionary picture 
and in a sequence where the strongest star formation occurs before strongest AGN activity, 
galaxies should follow a loop path moving up in the L^^? vs. Lagn diagram with the rise 
of star formation and t hen to the rig ht and perhaps down towards the most intense AGN 
phase (see discussion in iNetzerl ( 120091 )). We have placed SMGs with weak X-ray AGN from 
Alexander et al.l ( l2005al ) in Fig. [6] to show that their location in comparison to QSOs is 
consistent with such a path. SMGs without detected AGN would be found yet further to 
the left in Fig. [6l Comparing the location of individual submm-selected SMGs hosting weak 
AGN to the X-ray selected CDFS stacks has to be done with caution since the CDFS stacks 
themselves will contain similar SMGs contributing to the upturn in mm emission at high 
X-ray luminosity. A thorough analysis of possible evolutionary loops in such a L^i? vs. Lagn 
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diagram will thus have to be done on the basis of future individual measurements for sources 
in all the relevant parts of the diagram. 

4. Conclusions 

We have used the combination of the LESS 870//m survey and deep X-ray surveys of 
the (E)CDFS region to study star formation in the hosts of AGN covering a wide range of 
redshifts and luminosities centered around z~l and L2-wkeV ~ lO^^ergs"^. Stacking LESS 

data at the positions of all 895 AGN we detect at high significance submm emission at 
S87o^m=0.49±0.04mJy corresponding to average star formation rates of about SOMqYI^^ 
Using the good statistics to break down the sample according to AGN properties, we find 
an increase with redshift and little change with AGN luminosity, except for indications for 
an upturn in host star formation at L2-iQkeV > lO'^'^ergs^^. The bulk of the X-ray AGN do 
not show changes with AGN obscuration as expected from the merger evolutionary scenario, 
but such a behaviour may emerge at I/2-iofcey ^ lO^^ergs"^. 

Combined with results for higher luminosity AGN not properly sampled in the 0.25 
square degree ECDFS, we conclude that the bulk of deep survey X-ray AGN seem to be 
hosted by galaxies evolving secularly, with star formation rates similar to comparably massive 
non-active galaxies and no close link between AGN and global host star formation. In 
contrast, the most luminous L2-wkeV > lO^^ergs"^ AGN seem to follow a path where AGN 
activity and obscuration appear to be more closely linked to host star formation, likely via 
merger evolution. 

The properties of local Swift-BAT selected AGN with IRAS-hased far-infrared star 
forming luminosities are consistent with these two paths. The host star formation rates of 
moderate luminosity AGN are decreasing from z~l to z=0 similar to the decrease of SFR 
in non- active massive galaxies over this redshift interval. 
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and Sylvain Veilleux for discussions. D.L. thanks the Aspen Physics Center for hospitality 
during part of the preparation of this paper. I.R.S., K.E.K.C. R.J.I and J.W. acknowledge 
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Table 1. Stacking results for combined CDFS+ECDFS list 



Group -Nsource ^870fj,in (^map ^ subsample ^870t^m,residual LiR 

mJy mJy mJy mJy IO^^Lq 

(1) (2) (3) (4) (5) (6) (7) 



No grouping 

All CDFS & ECDFS 895 0.490 0.044 0.054 0.344 



Grouping by redshift 



z<1.2 


383 


0.279 


0.066 


0.083 


0.218 


1.15 


z>1.2 


365 


0.637 


0.068 


0.085 


0.445 


3.97 


z>2 


158 


0.597 


0.106 


0.129 


0.430 


3.79 


All z 


748 


0.453 


0.047 


0.059 


0.328 


2.48 


No good phot/spec z 


147 


0.695 


0.108 


0.134 


0.434 





Note. — Col. (1) — Group of X-ray AGN entering the stack. 
Col. (2) — Number of sources in stack. 

Col. (3) — Weighted mean of 870/im flux densities for stack. Measurements for 
stacks which are at < 3a, hence not significantly detected, are highlighted by italics. 
Col. (4) — Standard deviation of the weighted means, obtained by drawing many 
samples of ^ source random positions from the LABOCA map. 
Col. (5) — Standard deviation of the weighted means from drawing many subsamples 
of l^source from the full N=895 AGN sample. 

Col. (6) — Weighted mean of 870//m flux densities for stack on the residual map, 
obtained by subtracting all individual detections from the original map. For the stacks 
from the residual map, the corresponding errors a^ap and (Tsubsampie (not shown) are 
consistent with each other and ~12% less than amap for the full map. 
Col. (7) — Stacked IR luminosity, from the weighted mean of the luminosities of rest 
frame T=35K ^=1.5 greybodies individually matched to the redshift and 870//m flux 
of each X-ray AGN. Available only for samples with redshifts for each object. 
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Table 2. Stacking results for Tozzi et al. CDFS list 



Group -N^Source ^870iJ,m ^map ^subsample ^870iJ,m,residual L 



IR 



mJy mJy mJy niJy lO^^Lo 

(1) (2) (3) (4) (5) (6) (7) 



No grouping 



All CDFS 


302 


0.437 


0.074 0.094 


u.ozo 




Grouping by redshift 












Z<1 


145 


0.362 


0.109 0.134 


0.258 


1.47 


Z>1 


157 


0.506 


0.104 0.128 


0.394 


3.13 


z>2 


63 


0.779 


0.162 0.202 


0.561 


4.93 


Grouping by intrinsic 


rest frame 2-lOkeV luminosity 






L< 10^2grg g-i 


48 


0.379 


0.189 0.233 


0.289 


1.05 


L< lO^^erg 


138 


0.361 


0.111 0.138 


0.305 


1.52 


L> lO^^erg 


164 


0.500 


0.102 0.126 


0.348 


3.01 


L< lO^^erg 


260 


0.341 


0.086 0.101 


0.266 


1.71 


L> lO^^erg s"^ 


42 


1.036 


0.200 0.250 


0.718 


6.24 


Grouping by X-ray obscuring column from spectral fit 






Nh < lO^^cm-^ 


114 


0.315 


0.122 0.152 


0.280 


1.32 


Nh > 1022cm-2 


188 


0.511 


0.095 0.120 


0.357 


2.94 


Nh < 3 X 1022cm-2 


169 


0.421 


0.098 0.115 


0.325 


2.04 


Nh > 322cm-2 


133 


0.458 


0.113 0.140 


0.333 


2.66 


Nh < 1023cm-2 


225 


0.416 


0.087 0.108 


0.310 


2.08 


Nh > 1023cm-2 


77 


0.498 


0.149 0.123 


0.383 


3.01 


Sources with L2-10A;eV 


> 3 


lO^^ergs"^ 


^, grouping by X- 






ray column 












Nh < lO^^cm-^ 


18 


0.072 


0.316 0.378 


0.185 


0.40 


Nh > lO^^cm-^ 


80 


0.641 


0.147 0.180 


0.427 


3.95 


Nh < lO^^cm-^ 


48 


0.400 


0.189 0.233 


0.341 


2.42 


Nh > 102=^cm-2 


50 


0.667 


0.184 0.226 


0.422 


4.12 


Sources with L2-iokeV 


> 1 lO^^ergs-^ 


and spectroscopic 






redshift, grouping by X-ray column 








Nh < 1022cm-2 


6 


0.433 


0.531 0.653 


0.549 


2.23 


Nh > 1022cm-2 


21 


1.390 


0.285 0.352 


0.858 


8.43 


Nh < lO^^cm-^ 


13 


1.098 


0.366 0.447 


0.860 


6.40 
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Table 2 — Continued 



Group Ngource ^SYO/im ^map subsample ^870iim,residual L 



IR 



mJy mJy mJy mJy lO^^L© 

(2) (3) (4) (5) (6) (7) 



> 1023cm-2 14 1.244 0.346 0.432 0.719 7.57 

Sources with L2-iokeV < 3 x lO'^'^ergs"^ and spectro- 
scopic redshift, grouping by X-ray column 

Nh < 1022cm-2 96 0.360 0.133 0.166 0.298 1.47 

> 1022cm-2 108 0.415 0.124 0.156 0.307 2.23 

Sources with Lo.5-2fcey or L2-iokeV > 1 lO^'^ergs"^ and 

spectroscopic redshift, grouping by optical spectral type 

Type 1 (BLAGN) 10 0.737 0.410 0.508 0.743 4.29 

Type 2 17 1.436 0.315 0.388 0.815 8.64 



Note. — Definition of columns as in Table [H Stacks measured at < 3a, hence 
not significantly detected, are highlighted by italics. 
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Fig. 1. — Pixel flux histograms for the roughly 80% of the LABOCA 870yum map of the 
ECDFS with RMS less or equal 1.5 times the minimum noise of 1.07mJy. The black thick 
histogram represent s the full map, th e red thin histogram the residual map after subtraction 
of detected sources flWeiB et al.ll2009l ). Vertical lines indicate the mean flux of this region in 
the respective map. The blue thick histogram shows the flux distribution at the positions 
of the combined sample of 895 CDFS and ECDFS AGN. The blue thin histogram shows for 
easy comparison this AGN histogram scaled to the number of map pixels. 
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-50 50 -50 50 

Fig. 2. — Stamps showing the coadded LABOCA submm signal for all 895 CDFS and 
ECDFS X-ray sources in our combined hst. Left: Obtained using the total LABOCA flux 
map. The 27" FWHM beamsize of the LABOCA map is shown for comparison in the bottom 
left. Right: Obtained using the residual map after subtraction of all sources individually 
detected above 3.7a. The spatial scale is in arcseconds relative to the expected nominal 
center of the stacked beam. Contours are shown at 3, 6, 9 times the full stack noise of 
0.044mJy. The greyscale of the images runs from zero to the respective peak. 
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Fig. 3. — Top panel: X-ray selected AGN from the CDFS and ECDFS show an increasing 
trend of mean submm flux density with redshift. As in similar diagrams below, horizontal 
bars or arrows indicate the parameter range covered by a subsample, while vertical error bars 
indicate the la errors of the mean flux for this subsample. The symbols are placed at the 
median of the quantity on the abscissa for the respective bin - here redshift. The asterisk 
without redshift range indicates the mean flux of the entire combined sample. Results 
are shown separately for CDFS AGN only, and for the combined sample of CDFS and 
ECDFS AGN. Results for ECDFS sources without reliable redshift are shown at an arbitrary 
z, most likely these contain a significant fraction of presently unidentified high redshift 
sources. The m ean submm f luxes of (predominantly non-AGN) K<20 galaxies are added 
for comparison (IGreve et al.l |2009| . also based on LABOCA data). We also illustrate the 
submm flux expected from a greybody of given luminosity at different redshifts (dashed 
curve). Bottom panel: Results from stacking the IR luminosities of greybodies fitted to the 
redshift and map flux of each source. 
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Fig. 4. — Top panel: Mean submm flux density as a function of intrinsic rest frame 2-lOkeV 
luminosity of the AGN. The asterisk indicates the mean submm flux of the entire combined 
sample. The mean submm flux rises clearly above L2-iokeV = lO'^^ergs^^. These highest 
luminosities are also the sole driver of the flux increase in the wider L2-iokeV = 10^^ erg s~^ 
bin. Bottom panel: Results from stacking the IR luminosities of greybodies fitted to the 
redshift and map flux of each source. 
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Fig. 5. — Trends of mean submm flux density with AGN X-ray obscuring column density. 
Top panel: Sample split at A'^^^ = lO^^cm"^, as indicated by the dotted line. Bottom panel: 
Sample split at Nfj = lO^^cm"^. The asterisk shows the mean flux of the entire combined 
sample. We show results for the full CDFS sample, with no indication for a trend with 
Nh, and for the the most X-ray luminous CDFS sources with good redshifts. From the top 
panel, these may be more submm luminous above Nh ~ ^ "^^^ 
the difference is not significant in our sample. 



10 cm compared to below, but 
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Fig. 6. — AGN luminosities and star forming far- infrared luminosities for high redshift and 
local AGN. Moderate luminosity AGN have host star formation rather independent of AGN 
luminosity, but at a level rising with redshift. For high AGN luminosity, AGN luminosity and 
host star formation correlate. Red asterisks represent the four CDFS AGN stacks grouped by 
intrinsic 2-lOkeV X-ray luminosity, as also shown in FigJUand Table O The four stacks are 
for L2-iokeV < 10^^, < 10^^, > 10^^ and > lO^'^ergs"^ (see also the horizontal bars indicating 
the range of AGN luminosities entering) and are plotted at the median AGN luminosity 
of the sources in each stack. The diagonal dotted green line indicates the relation for a 
continuous host and black hole growth, reflecting the local universe relation between black 
hole mass and bulge mass. The black square represen ts the comb i ned X -ray obscured and 



unobscured samples of more luminous X-ray AGN from iPage et al.l (120041 ) and IStevens et al. 
(120051 ) [P04, SOS], see also Section 3. AGN and star formation luminosities have been 
obtained analogously to the CDFS stacks. The pink ellipse indicates the location in such a 
diagram of mm-bright and optically very luminous high-z QSOs studied by lLutz et al.l (120081 ) 
[L08]. The pink diamond reflects the suggested median location of the parent pop u lation 
of these luminous high-z Q SOs including mm- fai nt o nes, as discuss e d in Lutz et al.l (120081 ) 



and based on the surveys of IPriddev et all (l2003f ) and lOmont et all (l2003f ) [P03, O03]. The 



dashed line is a simple constant plus power law relation between AGN and star forming 
luminosity approximating the high redshift AGN. We also in clude local un i verse (z< 0.3) 
AGN, the blue ellipse representing local PG QSOs studied by iNetzer et al.l (120071 ) and the 
blue asterisks stacks of local hard X-ray selec ted AGN based on Swift-BAT and IRAS data 
(Section 3.3). The blue square represents the lElvis et al.l ( 1l994j ) QSO sample at its median 
AGN luminosity converted to our cosmology, with the error bar indicating the la variation 
in FIR/optical (star formation/ AGN) flux ratio of obj ects in that sample. Submillimeter 
galaxies with weak X-ray AGN (lAlexander et al.ll2005al ) are shown for comparison. 



